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ABSTRACT 



Dynamic instability of planing craft on calm water, porpoising, 
is a phenomenon which has not been properly understood. Empirical 
relations are available for predicting the regime of stability. The re- 
lations, when compared, lead to conflicting design requirements to in- 
crease stability. 

It is therefore desirable to develope a theoretical approach to 
the problem so that the effects of beam, deadrise angle, etc. on 
stability can be studied. 

The results of the investigation imply that a decrease in dead- 
rise angle, a decrease in beam and an increase in distance from LCG 
to transom result in an increase in stability. Changes in shaft angle 
and vertical height of the center of gravity and moment of inertia have 
very little effect on the stability of a boat while it is planing. However 
further investigation is required to verify these results. 

In conjunction with this paper, a computer program was written 
which can be used in the design of planing craft to predict boat attitude, 
wetted surface area, drag and effective hor3e power. This program 
will be available for use in the XIII Department library. 
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I. INTRODUCTION 

The unstable motions of planing craft have been under study for 
many years and have been the subject of much literature (see biblio- 
graphy). The ability to be able to predict the stability characteristics 
of a particular hull in the early stages of design is of importance to 
naval architects. A knowledge of the effects of variables such as beam, 
deadrise angle, etc. on stability would permit intelligent corrective 
action to be taken to increase the dynamic stability of existing craft. 

The problem of planing craft stability involves many variables 
and empirical relations between some of the design variables have 
been developed to predict dynamic stability. 

Two formulae recently developed emperically from expfrimental 
data, (2) and (12), result in conflicting design requirements to increase 
stability (see Appendix C). It is therefore desireable to develope a 
theoretical approach to the problem so that the effects of design var- 
iables can be determined independently of experimental data. 

Perring (10) attempted a theoretical approach. His lack of 
success can be attributed to a number of causes. The foremost of 
these being lack of sufficient experimental and theoretical information 
to predict the stability derivatives accurately and the ommission of 
important terms. 
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II. THEORY 

A planing hull, as a rigid body, has six degrees of freedom. This 
study treats the boat as a two degree of freedom system by investigat- 
ing what is considered to be the most important motions, heave in the 
Z -direction and pitch about the Y-axis. 

The equations of motion are nonlinear. To facilitate the solu- 
tion of these equations it is necessary to linearize them. 

Linearized equations of motion for ships have been developed 
by Abkowitz (1), Korvin-Korvosky (7) and others. Those of Abkowitz 
are most complete. If the coefficients, i. e. stability derivatives, are 
substituted into these equations and then the result transformed into 
the frequency domain, it should be possible to evaluate the stability 
by using the Routh criterion. ( 5 ). 

The method used to predict stability or lack thereof proceeds 
as follows: 

1. The stability derivatives are determined, see Appendix A. 

2. The stability derivatives are substituted into the linearized equations 
for ship motion. (1) 

3. The resulting equations are transformed by substitutions of the form 
z * Zmax e and 0 * 8max e . 

4. An equation in S is obtained. 

5. The Routh discriminant is evaluated for the fourth order equation 
in S, see Appendix A. 
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III. DESCRIPTION OF WORK ACCOMPLISHED 

The firststep toward a solution was to determine the stability 
derivatives and combine them to form the coefficients of the linearized 
equations of motion, see Appendix A. The coefficients were then non- 
dimensionaiized using beam as the non-dimenstonalizing length (12). 

A computer program was written to solve for the Routh discrim- 
inants, see Appendix H, using as input data the results from a series 
of tests run at the David Taylor Model Basin (2), see Table 3. The 
resulting discriminants were then plotted against speed showing a 
consistency in the directions of the paths, however there was no ob- 
vious difference between stable and unstable* boats, see Figure 9. 

At this point, an attempt was made to determine the roots of the 
fourth order stability equation to examine their loci using a computer 
program from the MIT "SHARE" library. (SHARE No. 1514 RTSCH). 
RTSCH proved to be unsatisfactory. The answers obtained from this 
program are seriously in error even for the simplest of input equa- 
tions . 

It was then decided to vary in turn what seemed to be the most 
important variables: DIFB1, DIFC1, B2, Dl, D2, El, Gl, and G2. 

This was done to determine the effect of changes in their magnitude 
on the Routh discriminant. From Figure 10 it can be seen that vary- 
ing G2 roughly grouped the stable and unstable boats with the unstable 
group centered about G2 equal to 0. 38 x G2 at the point of zero Routh 
discriminant. The program was then run with G2 equal to 0. 38 x G2 
so that the loci of the discriminants could be examined. The results 
are shown in Figure 11. Based on these results it was decided to in- 

1. Unstable boats, as refered to in this paper, are those which porpoised 
at a Fv less than 6. 0. Stable boats are those which had not porpoised 
before maximum test speed (2) was attained (Fy * 6. 0). 
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vestigate G2 further. Since the term G2 is made up of two parts, force 

£{arm) , a(force) .. . . . . . x 

x -r + arm x — - , it was decided to vary these two parts 

cZ 02 ^ 

independently to see if better correlation could be achieved in either 
of the two groups. Correlation was not improved, see Figures 12 
through 17. However it was observed that the discriminants for the 
stable boats changed very little with changes in G2, or its parts. On 
the other hand the discriminants of the unstable boats changed a great 
deal with changes in G2. This lead to the conclusion that there must be 
a term in the coefficients of the characteristic equation which over- 
powered G2 when the boat was stable but was of the same magnitude 
as G2 when the boat was unstable. 

Based on information obtained thus far, it was decided to investigate 
each term of the Routh discriminant to see which of the coefficients were 
controlling for the stable and unstable boats. Values of the coefficients 
obtained from program 1 in Appendix H were inserted into each term 
manually and inspection of the results was unfruitful. No obvious dif- 
ference could be detected between stable and unstable boats. It was 
concluded that the interaction was much more subtle. 

A closer inspection of each term indicated that the coefficients 
and , D1 and D2, may interact with important effects and this became 
the final step in the investigation of the coefficients. The results are 
shown in Figures 18 through 20. At this point the investigation of the 
coefficients of the equations of motion was terminated because of time 



1. At this point it was necessary to reduce the number of plots to 
three stable and three unstable boats in order both to simplify the plots 
and to make more efficient use of computer time. The unstable boats 
were selected by choosing two which had axis intercepts fairly close 
together and third whose intercept was remote from these (models 
4665-3, 4666-13, and 4668-9 in Figure 10). 



